Phase-space reconstruction of an atòmic chaotic system 
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We consider the dynamics of a single atom submitted to periòdic pulses of a far-detuned standing 
wave generated by a high-finesse optical cavity, which is an atòmic version of the well-known "kicked 
rotor". We show that the classical phase-space map can be "reconstructed" by monitoring the 
transmission of the cavity. We also studied the effect of spontaneous emission on the reconstruction, 
and put limits to the maximum acceptable spontaneous emission rate. 

PACS numbers: 42.65.Sf, 03.75.Dg, 39.20.+q 32.80.Pj 



I. INTRODUCTION 

The kicked rotor (KR) is a widely explored system that 
has a paradigmàtic status in studies of classical [| and 
quantum chaos In recent years, the advent of làser 

cooling of àtoms allowed the realization of an atòmic ver- 
sion of the KR HM, which has been used in experiments 
by many groups 0,0, 0, @ • This consists simply in plac- 
ing laser-cooled àtoms in a pulsed làser standing wave. 
In adequate (but rather general) conditions, the atòmic 
center of mass motion "feels" the presence of the standing 
wave as an "optical potential" which varies sinusoidally 
in space. One can then easily modulate (with period T) 
the light intensity in the form of short pulses to obtain a 
hamiltonian of the form 
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n=0 



where M and p are resp. the atom mass and center 
of mass momentum, Vq is the amplitude of the optical 
potential (which is proportional to I/A, where / is the 
intensity of the standing wave and A the laser-atom de- 
tuning), fci = 2-ïï/Xl the wave number of the standing 
wave (in the x— direction), and S T a square function of 
duration r and amplitude 1. In the limit t — > 0, or, more 
exactly, in the limit 
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one can put S t /t — > S and one obtains exactly the 
hamiltonian of the KR. In order to obtain this hamilto- 
nian in its Standard form, we perform the normalizations 
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where the only free parameter is the normalized kick 
strength K = W k 2 L TT/M. 

Many of the recent experimental studies were aimed 
at quantum aspects of the KR dynamics, as "dynamical 
localization" 0, EJ , quantum resonances , or "chaos- 
assisted" tunneling [D], H2 • Quantum effects in such sys- 
tem manifest for times larger than the "Ehrenfest time", 
when quantum dynamics begins to take over classical dy- 
namics. This is the domain of "quantum chaos", which is 
defined as the quantum dynamics of systems which are 
chaotic in the classical limit. It is well known that the 
quantum behavior of such systems bears little relation to 
the classical chaotic behavior of their classical counter- 
parts, basically because Schròdinger equation is linear. 
One of the most exciting qüestions about quantum chaos 
is "How the initially classical dynamics and latter quan- 
tum dynamics diverge from one another?" Clearly, the 
critical point in this divergence is situated between the 
Ehrenfest time tg ~ Tln(l/£) and the so-called "localiza- 
tion time " tj ~ T(K/k) 2 , where k is the effective Planck 
constant [19| . Fortunately, in the atòmic KR the param- 
eters can be widely changed which allows one to adjust 
these times to experimentally convenient seales. How- 
ever, even if the adequate time domain is accessible, it is 
still dimcult to put into evidence the transition between 
classical and quantum regimes because, usually, the kind 
of measurements performed on a microscopic system is 
different from - and often not simply related to - the 
kind of measurement performed on macroscopic systems. 
Typically, studies of classical dynamics rely on the no- 
tion of trajectory (phase-space maps, first-return maps, 
Poincaré sections...) whereas studies of quantum dynam- 
ics emphasize mean quantities and probability distribu- 
tions. The purposc of the present paper is to describe and 
study an experimentally realizable situation in which re- 
construction of the atom center of mass dynamics can 
be done, thus reconciliating quantum- and classical-type 
measurements. Hopefully, such a setup would provide a 



better understanding of the above-mentioned transition 
between classical and quantum dynamics. Throughout 
this paper, we shall treat the atòmic motion as a classi- 
cal motion: we shall consider that the atòmic center of 
mass position and momentum are enough to character- 
ize the dynamics. This is true as long as the dynamics 
is essentially classical, that is, for times shorter than the 
localization time, or as far as the sources of decoherence 
- mainly spontaneous emission in the present case - act 
often enough to provoke a quantum coherence col·lapse 
before interference effects can develop and appreciably 
change the classical behavior. 



II. KICKED ROTOR DYNAMICS 

In recent years, the quality of optical cavities has been 
dramatically improved. It is now possible to build cav- 
ities of very high finesse [ÏÜL lïij ] . In such a cavity, the 
index of refraction generated by one atom is enough to 
change mesurably the resonance condition, and the cav- 
ity can thus be used to "count àtoms one by one" [Ï5| . 
As an atom moves inside the cavity, it "probes" regions 
where the radiation intensity is not the same. For exam- 
ple, cavity modes usually have a gaussian profile, whose 
intensity is larger at the center than in the border. In the 
same way, as the field inside the cavity forms a standing 
wave, changing the longitudinal position of the atom also 
changes the intensity to which it is submitted (Fig. [TJ. 
This is the situation we shall consider in the present pa- 
per. Because of the Kerr effect, the refractive index of the 
atom is modified by the intensity of the radiation to which 
it is exposed, and the changing of the refractive index in 
turn changes the cavity resonance condition. By measur- 
ing the cavity transmission one can extract information 
on the position of the atom. This is an "autoconsisteney" 
problem: the change in the resonance condition changes 
the cavity field, that in turn changes the atom refractive 
index. This means that the amplitude of the kicks felt 
by the atom will also depend on its position, which is in 
contradiction with the Standard definition of the KR. We 
shall take this fact into account in our simulations, but 
will arrange parameters in order to make the variation 
of the kick strength K small enough to not disturb too 
much the dynamics of the KR. 

An experimental limitation that must also be taken 
into account is that one cannot change instantaneously 
the radiation intensity inside a cavity. The "lifetime" of a 
photon inside the cavity is roughly r c — L/(cT), where T 
is the (intensity) transmission cocfficicnt of the mirrors 
and L the length of the cavity. The duration r of the 
kicks shall thus be at least of a few r c . However, the 
validity condition for identifying the kicks with a delta 
function, Eq. @, must also be satisfied. This implies 
that 




Figure 1: Schematic view of the proposed setup. If the atom 
is at a crest (a) or at a node (b) of the cavity electromagnètic 
field, the Kerr nonlinear refractive index produced is differ- 
ent, the resonance condition of the cavity is modified, and 
the transmitted intensity changes, allowing detection of the 
longitudinal position of the atom. 



1/2 

where we used the order of magnitude value (p 2 ) <~ 

2hkh\TN after TV kicks are performed |2íj. Thus given 
T, there is a maximum number of kicks (i.e. a maximum 

duration of the experiment), N max ~ \\ 2 L M cT / {2hL)\ . 
For cesium this gives N max ~ (T/L) 2 x 10 9 m 2 . This 
condition is not hard to satisfy: for T = 1CP 6 and L = 
1CP 4 m, N max ~ 10 5 , whereas experiments currently last 
only for a few hundred kicks. 

Integrating the equations of motion obtained from 
Eq. ® over oneperiod T = 1 produces the so-called 
"standard map" [lj 



X n +i — X n + P n 



Pn+l — Pn — K n+ i sin X n+ i. 



(5a) 



(5b) 



where we took into account the fact that, due to the 
nonlinear Kerr effect, K n = K(X n ): the intensity of the 
kick depends on the position of the atom. In constructing 
the phase-space portrait of the kicked rotor, it is usual 
to simplify the display by taking all the quantities "mod- 
ulus 27r" (we shall represent the modulus 27r of a given 
quantity A by A[27r]). This is due to the fact taht the 
potential is periòdic in space, of period 27r in reduccd 
units. This means that if P n — l-nq + p, with p < 2ir and 
q integer, then, after Eq. (JïïaJ, X n+ i — X n + p + 2nq. 
Thus A„ + i[27r] = A„[27r] + p: the physics is the same 
for momenta differing by an integer múltiple of 27r. We 
shall thus, according to the usual convention, draw phase- 
space maps by plotting A[27r] vs. P[2n]. We show a 
typical plot of the KR phase-space in Fig. 0Ja). 



III. CAVITY DYNAMICS 
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A schematic view of the proposed experiment is shown 
in Fig.^ The atom is placed inside a high-fmesse cavity, 



3 



and its presence shifts the resonance condition of the cav- 
ity in a nonlinear way: this shift depends on the radiation 
intensity seen by the atom, and thus on its longitudinal 
position. The intensity transmitted through the cavity 
can thus be related to the position of the atom. However, 
this position is defined only with respect to the closest 
node; shifting the position of the atom by a period of 
the standing wave produces the same transmission signal 

m 

If a làser beam of wavenumber L·l intensity Iq is in- 
jected into an empty Fabry-Perot cavity made of two 
identical mirrors a distance L/2 apart, of transmission 
coeficient T and reflection coeficient 1Z — 1 — T, the in- 
tracavity intensity is given by the well-known Airy func- 
tion 

l/T 



L = 



l + i 72 sin 2 (^J^ 

leading to a resonance Rabi frequeney 
fi 2 = ÍIq cos 2 (L·lx) 

( 37rr \ l/T 
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where F = 2\/lZ/T is the finesse of the cavity divided 
by 7r/2 and T the natural width of the atòmic transition. 
The parameter K [Eq. J3] turns out to be 0: 
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where I s = hckfT/(6ir) is the saturation intensity ~ 2 
mW/cm 2 for cesium, uj r = hk L /(2M) is the recoil fre- 
queney (u) r /(2ir) ~ 2 kHz for cesium) and A the laser- 
atom detuning. 

In Eq. JHJ), L op t = L + AL is the optical length of the 
cavity, L is its physical length and AL is the correction 
due to the presence of the atom: 
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where L at is the linear size of the effective volume oceu- 
pied by the atom. The definition of this volume depends 
on the details of the experiment being performed. In the 
case of laser-cooled àtoms, we can take this volume as 
the volume of the atom cloud in a magneto-optical trap, 
roughly L at ~ 100 /zm, which gives Ziï/[2(kL Lat) 2 ] ~ 
10~ 5 . For A rí 10 3 r, which is a typical value, the de- 
phasing kLAL induced by a lone atom is roughly 10~ 8 , 
which is hard to detect directly. The effect is however 
enhanced by the presence of the cavity. 

The detected signal is the intensity transmitted 
through the cavity, given by 

1 - . , ■ h l/T @y 



1 + F 2 sin" 
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which depends on the radiation intensity seen by the 
atom through L opt , and thus on the atom position. If 
F 3> 1, the half-width of the Airy peak corresponds 
roughly to 



F 2 sin 2 
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One easily finds 
AI t 



^± \k L AL\ = | \k L AL\ 



where we supposed fc^L w 2/F . Then, taking kLAL « 
10 -8 , an experimentally realistic value 0] of F ~ 10 6 
produces a detectable 1%-variation of the transmitted 
intensity p^. 

Techniques based on the Kerr nonlinearity similar to 
that proposed here have been recently used for monitor- 
ing the radial motion of an atom trapped in high-finesse 
cavity 11311141 . In our case, it is the longitudinal motion of 
the atom that should be monitored. This is complicated 
by the periòdic character of the standing wave inside the 
cavity. We shall discuss in the next section the algorithm 
we developed to deal with this problem. 



IV. THE RECONSTRUCTION ALGORITHM 

In order to reconstruct the phase-space map from the 
cavity transmission signal, we suppose that the trans- 
mission value corresponding to each kick n is recorded, 
forming a temporal series S n . The information given by 
the bare signal is however incomplete. Fig.|2]displays the 
first return map corresponding to the temporal series S n . 
Its topology do not evoke at all the KR phase-space map 
shown in Fig. Let us consider the information that 
can be extracted from two successive vàlues of the cavity 
transmission, Si at kick 1 and S2 at kick 2. Referring 
to Fig. |3 we see that the cavity transmission leaves an 
indeterminacy on the position, as the points X\ , X[ and 
X'/ correspond to the same transmission. Generally two 
classes of points lead to the same transmission 



X and X' = X + 2ttç 

X and X" = 2nq - X. 
(q is an integer). 
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Figure 2: First return map S„ x S n +i for the cavity- 
transmission signal, with í 1 = 10 6 . Structures are clearly 
visible, but no direct relation to the topology of the phase- 
space of the KR (Fig. can be drawn. 
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Figure 3: Correspondence between the transmission signal 
and the atom longitudinal position inside the cavity. 



X Xi X\ + 2n qi 2nq! - Xi 



x 2 


X2 — X\ 


X 2 - Xi - 2ir qi 


X 2 + Xi - 27rgi 


X 2 + 27TÇ2 


X 2 - Xi + 2nq 2 


X2-X 1 + 2TT(q 2 -qi) 


X 2 +Xi +27r(q 2 - gi) 


2-KQ2 — X2 


-(X 2 +Xl) +27VQ2 


-(X 2 +X 1 ) + 2iv(q2-qi) 


-(X 2 - Xi) + 27r(g 2 - gi) 



Table I: Possible momentum (X 2 — Xi) determinations for different choices of the position determination (çi and g 2 are 
integers). Eliminating the terms of the form 2-7T x (integer number), the only possibilities are X2 + xi, — (22 + Xi), %2 — xi 
and — (x2 — xi), with Xi — Xj[27r]. Considering that trajectòries corresponding to momenta P and —P are identical in the 
phase-space map leaves only two choices, X2 + x\ and X2 —x\. 



However, as it is X[27r] that is displayed in the phase- 
space map, the value of the integer number q in Eqs. H12(l 
is irrelevant. This leaves us with just two possible vàlues, 
x and 2n — x, where x = X[2tt]. As the motion of the 
atom is frec bctwccn two kicks, the momentum associated 
to the kick n is P n = X n+ i — X n . Moreover, we note 
that P and — P lead to the same trajectory in the phase 
portrait. Introducing the relevant (for phase-portrait- 
plotting purposes) value p n = P n [2ir] , wc find that there 
are only two distinct choices for the momentum: pi — 
x 2 — X\ or pi — x 2 + X\ (see table[IJ|. We must thus chose 



arbitrarily one of the two. If we now include the next 
value X3 of the position (with its own arbitrary choice 
of the determination), we can deduce in the same way 
a value for p 2 . But Eq. (|5b|) implies that the successive 
vàlues p\ and p 2 of the momentum must satisfy 

P2 =Px — {K 2 siax2) [2tt], (13) 

which allows us to test the coherence of the choices. 

In practice, our algorithm proceeds as follows: i) For 
each point S n in the time series of cavity transmission, we 
first estimate a value for the intracavity intensity 1^ us- 
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Figure 4: Comparison between the original (a) and the recon- 
structed (b) KR phase space-map. Parameters are K — 0.8, 
F = 10 6 . Larger structures are easier to reconstruct than 
smaller ones. 



ing Eq. (JüJ with L opt — L. This value of Ic°^is then used 

to estimate the atom position Xn^ and to calculate a first 
correction to L, = L + ALi, which in turn is used 
to determine a corrected value Ic and a corrected value 
of the position Xn ■ The process is iterated until it con- 
verges to produce the value of the position correspond- 
ing to the n th kick, x n . Ten steps are usually enough to 
insure the convergence of the process. ii) We choose ar- 
bitrarily between the two possible determinations x n or 
2tï — x n of the position and calculate the corresponding 
value of p n (choosing arbitrarily between x n — x n -\ and 
x n + x n -i). iii) Each time three successive momentum 
vàlues havc bccn dctcrmincd, we use Eq. (JT3J to test the 
coherence of the choices. If this condition is violated, we 
recalculate the last three points with a different choicc 
of the position determinations. iv) Once a trajectory is 
reconstructed, one starts again with another time series 
corresponding to another set of initial conditions. 

Fig. 21 compares the exact KR phase-space map with 
the reconstructed one. We see that the agreement is 
good. The number of restarts due to the coherence test 
is 380 over 3000 data points. 



V. EFFECT OF THE SPONTANEOUS 
EMISSION 

The good reconstruction obtained in the preceding sec- 
tion is however impossible to achievc in a real experiment, 
because of the unavoidable presence of spontaneous emis- 
sion in the system. However, the spontaneous emission 
rate can be controlled - to a certain extent - in such sys- 
tems, because it roughly seales with the inverse of the 
square of the laser-atom detuning, whereas the optical 
potential amplitude Vq seales with the inverse of the de- 
tuning; one can thus adjust the parameter in such way 
that spontaneous emission level is tolerable during the 
experiment duration. In the one dimensional situation 
we are considering, the effect of spontaneous emission is 
to add to the value of p in the average a unit of photon 



momentum zt^fc^ (corresponding to k/2 = 2hkf j T/M in 
reduced units), with an arbitrary sign, each time it hap- 
pens. We introduced spontaneous emission in our simu- 
lations by a simplified Monte-Carlo procedure. For each 
pulse of duration r, we calculate the probability for a 
fluorescence cycle to happen 
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We suppose that this probability is less than one, that 
is, there is at most one fluorescence cycle per pulse. A 
random number < a a < 1 is picked, and if a a < g a 
the atom performs a fluorescence cycle. In such case, 
another number rj a taking randomly one of two possible 
vàlues —k/2 and +k/2 is picked to decide of the direc- 
tion of the recoil of the atom [23j . The atom decay back 
to the ground state produces a new recoil in a random 
direction given by another random number r] sp = ±k/2. 
We suppose that T 3> 24], so that, with respect 

to the atom dynamics, the whole fluorescence cycle can 
be considered as instant aneous. The effect of the spon- 
taneous emission during the n th kick is thus to change 
the momentum p n — > p n + (r] a +r) sp ). This change in 
the momentum displaces the point corresponding to the 
p n in the phase-space of a quantity that can be as large 
as £[27r], compared to the width of the map 2tt. The 
effect of spontaneous emission thus increases as the ef- 
fective Planck constant k of the system increases, that 
is, as the system becomes more "quantum". Fortunately, 
the classical type of dynamics we are interested in here 
corresponds to fc < 1. 

Fig. [S] compares the phase-space reconstruction ob- 
tained with different vàlues of g a . The number of restarts 
due to bad choices of position/momentum determina- 
tions was typically two times larger than in the absence 
of spontaneous emission. We can estimate the maximum 
spontaneous emission rate per kick still allowing the re- 
construction of the phase map as g ma x ~ 0.05, which 
means that, in order to perform this experiment, the pa- 
rameters should be chosen so that the localization time 
be of the order of Íl ~ 9 ma x ~ 20 kicks. However, 20 
kicks are not enough to produce a good phase space por- 
trait; a better combination of parameters would be, e.g., 
g ~ 0.01 and íl ~ 100. This vàlues are still compatible 
with current experimental vàlues. 



VI. CONCLUSION 

We have proposed and numerically tested a method 
allowing reconstruction of the classical dynamics of the 
"atom-optics" kicked rotor formed by an atom exposed 
to pulses of an intracavity standing wave. We have dis- 
cussed the relevant parameters of the system and showed 
that they can fit with real experimental situations. Our 
method can be used to probe interesting phenomena, like 
the transition from the classical dynamics to the quan- 
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Figure 5: Influence of the spontaneous emission in the phase- 
space reconstruction. Parameters are K = 0.8 and F — 10 6 . 



tum dynamics as the spontaneous emission rate becomes 
smallcr than the inverse of the localization time. 



It is rather difhcult to guess what will happen when 
this limit is approached. In the quantum case, the atom 
cannot be characterized by a single point in the phase 
space. This means that the spreading of the wavefunction 
must be considered, and the classical notion of trajectory 
is not any more a good one. In particular, the the cavity 
transmission will be affected by the spatial spreading, as 
it allows the atom to "probe" different intensity regions. 
One must also consider the decoherence effect of spon- 
taneous emission that reduces the wavepacket if it hap- 
pens often enough. An interesting way to explore such 
situation is to compare a Wigner-function picture of the 
quantum phase-space in presence of decoherence with the 
classical map [|ï| . If the present method proves able to 
experimentally explore the classical/quantum transition, 
it would be highly valuable for studies of the nontrivial 
relation between classical and quantum chaos. 
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